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Solids produced as a result of a fast quench across a freezing 
or a structural transition get stuck in long-lived metastable 
configurations of distinct morphology, sensitively dependent 
on the processing history. Martensites are particularly well 
studied examples of nonequilibrium solid-solid transforma- 
tions. Since there are some excellent reviews on the subject, 
we shall, in this brief article, mainly present our viewpoint. 



I. NONEQUILIBRIUM STRUCTURES IN SOLIDS 

What determines the final microstructure of a solid 
under changes of temperature or pressure ? This is an 
extremely complex issue, since a rigid solid finds it diffi- 
cult to flow along its free energy landscape to settle into a 
unique equilibrium configuration. Solids often get stuck 
in long-lived metastable or jammed states because the 
energy barriers that need to be surmounted in order to 
get unstuck are much larger than fc^T. 

Such nonequilibrium solid structures may be obtained 
either by quenching from the liquid phase across a freez- 
ing transition (see Ref. [n| for a comprehensive review), 
or by cooling from the solid phase across a structural 
transition. Unlike the former, nonequilibrium structures 
resulting from structural transformations do not seem to 
have attracted much attention amongst physicists, apart 
from Refs. MM, possibly because the microstructures and 
mechanical properties obtained appear nongeneric and 
sensitively dependent on details of processing history. 

Metallurgical studies have however classified some of 
the more generic nonequilibrium microstructures ob- 
tained in solid (parent/ austenite) - solid (product/ fer- 
rite) transformations depending on the kind of shape 
change and the mobility of atoms. To cite a few : 

• Martensites are the result of solid state transfor- 
mations involving shear and no atomic transport. 
Martensites occur in a wide variety of alloys, poly- 
meric solids and ceramics, and exhibit very distinct 
plate-like structures built from twinned variants of 
the product. 

• Bainites are similar to martensites, but in addition 
possess a small concentration of impurities (e.g. 
carbon in iron) which diffuse and preferentially dis- 
solve in the parent phase. 



• Widmanstatten ferrites result from structural 
transformations involving shape changes and are 
accompanied by short range atomic diffusion. 

• Pearlites are a eutectic mixture of bec Fe and the 
carbide consisting of alternating stripes. 

• Amorphous alloys, a result of a fast quench, typi- 
cally possess some short range ordering of atoms. 

• Poly crystalline materials of the product phase are 
a result of a slower quench across a structural tran- 
sition and display macroscopic regions of ordered 
configurations of atoms separated by grain bound- 
aries. 

That the morphology of a solid depends on the de- 
tailed dynamics across a solid-solid transformation, has 
been recognised by metallurgists who routinely use time- 
temperature-transformation (TTT) diagrams to deter- 
mine heat treatment schedules. The TTT diagram is 
a family of curves parametrized by a fraction 5 of trans- 
formed product. Each curve is a plot of the time required 
to obtain 5 versus temperature of the quench (Fig. 1). 
The TTT curves for an alloy of fixed composition may 
be viewed as a 'kinetic phase diagram'. For example, 
starting from a hot alloy at t = equilibrated above the 
transition temperature (upper left corner) one could, de- 
pending on the quench rate (obtained from the slope of 
a line T(t)), avoid the nose of the curve and go directly 
into the martensitic region or obtain a mixture of ferrite 
and carbide when cooled slowly. 
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Fig. 1 TTT curves § for steel AISI 1090 (0.84% C + 0.60% 
Mn). A: austenite (fee), F: ferrite (bee), C: carbide (FesC). 
Curves correspond to 0, 50 and 100 % transformation. Below 
a temperature M s , the metastable martensite (M) is formed 
- the transformation curves for martensites are horizontal. 

It appears from these studies that several qualitative 
features of the kinetics and morphology of microstruc- 
tures are common to a wide variety of materials. This 
would suggest that there must be a set of general princi- 
ples underlying such nonequilibrium solid-solid transfor- 
mations. Since most of the microstructures exhibit fea- 
tures at length scales ranging from 100A-100/im, it seems 
reasonable to describe the phenomena at the mesoscopic 
scale, wherein the solid is treated as a continuum. Such 
a coarse-grained description would ignore atomic details 
and instead involve effective continuum theories based on 
symmetry principles, conservation laws and broken sym- 
metry. 

Let us state the general problem in its simplest con- 
text. Consider a solid state phase diagram exhibiting 
two different equilibrium crystalline phases separated by 
a first order boundary (Fig. 2). An adiabatically slow 
quench from Ti n — > Tfi„ across the phase boundary in 
which the cooling rate is so small that at any instant 
the solid is in equilibrium corresponding to the instanta- 
neous temperature would clearly result in an equilibrium 
final product at Tfi n . On the other hand, an instanta- 
neous quench would result in a metastable product bear- 
ing some specific relation to the parent phase. The task 
is to develop a nonequilibrium theory of solid state tran- 
formations which would relate the nature of the final ar- 
rested state and the dynamics leading to it to the type 
of structural change, the quench rate and the mobility of 
atoms. 

In this article we concentrate on the dynamical and 
structural features of a class of solid-solid transforma- 
tions called Martensites. Because of its commercial im- 
portance, martensitic transformations are a well studied 
field in metallurgy and material science. Several clas- 
sic review articles and books discuss various aspects of 
martensites in great detail p]. The growing literature 



on the subject is a clear indication that the dynamics of 
solid state transformations is still not well understood. 
We would like to take this opportunity to present, for 
discussion and criticism, our point of view on this very 
complex area of nonequilibrium physics fifffl . 

We next review the phenomenology of martensites and 
highlight generic features that need to be explained by a 
nonequilibrium theory of solid state transformations. 



II. PHENOMENOLOGY OF MARTENSITES 

One of the most studied alloys undergoing martensitic 
transformations is iron-carbon J5l . As the temperature is 
reduced, Fe with less than 0.02% C undergoes an equilib- 
rium structural transition (Fig. 2) from fee (austenite) to 
bec (ferrite) at T c = 910°C . An adiabatic cooling across 
T c nucleates a grain of the ferrite which grows isotrop- 
ically, leading to a polycrystalline bec solid. A faster 
quench from T;„ > T c to Tfi n < M s < T c (where M s : 
martensite start temperature) produces instead a rapidly 
transformed metastable phase called the martensite, pre- 
empting the formation of the equilibrium ferrite. It is 
believed that martensites form by a process of hetero- 
geneous nucleation. On nucleation, martensite 'plates' 
grow radially with a constant front velocity ~ 10 5 cm/s, 
comparable to the speed of sound. Since the transfor- 
mation is not accompanied by the diffusion of atoms, 
either in the parent or the product, it is called a dif- 
fusionless transformation. Electron microscopy reveals 
that each plate consists of an alternating array of twinned 
or slipped bec regions of size « 100A. Such martensites 
are called acicular martensites. 
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Fig. 2 Phase diagram of Fe-C (weight percent of C < 
0.02%). M s is the martensite start temperature. 

The plates grow to a size of approximately \\xm before 
they collide with other plates and stop. Most often the 
nucleation of plates is athermal ; the amount of marten- 
site nucleated at any temperature is independent of time. 
This implies that there is always some retained fee, parti- 
tioned by martensite plates. Optical micrographs reveal 
that the jammed plates lie along specific directions known 



as habit planes. Martensites, characterised by such a 
configuration of jammed plates, are long lived since the 
elastic energy barriers for reorganisation are much larger 
than ksT . 

A theoretical analysis of the dynamics of the marten- 
sitic transformation in Fe-C is complicated by the fact 
that the deformation is 3-dimensional (Bain strain) with 
3 twin variants of the bcc phase. Alloys like In-Tl, In-Pb, 
Mn-Fe and high-T c ceramics however, offer the simplest 
examples of martensitic transformations having only two 
twin variants. For instance, the high-T c cuprates un- 
dergo a tetragonal to orthorhombic transformation J2|. 
The orthorhombic phase can be obtained from the tetrag- 
onal phase by a two-dimensional deformation, essentially 
a square to rhombus transition. Experiments indicate 
that all along the kinetic pathway, the local configura- 
tions can be obtained from a two-dimensional deforma- 
tion of the tetragonal cell. This would imply that the 
movement of atoms is strongly anisotropic and confined 
to the ab-plane. Thus as far as the physics of this trans- 
formation is concerned, the ab-planes are in perfect reg- 
istry (no variation of the strain along the c-axis). In 
the next two sections we shall discuss our work on the 
dynamics of the square to a rhombus transformation in 
2-dimensions using a molecular dynamics simulation and 
a coarse-grained mode coupling theory. 



tial block of size £ and a time interval r, from the in- 
stantaneous positions u of the particles. These include 
the deformation tensor dui/dx k , the full non linear strain 
eij , and the vacancy field <j> = p — ~p (p = coarse grained 
local density, ~p = average density). We have kept track 
of time dependence of these coarse-grained fields during 
the MD simulation. 



LIQUID 




Fig. 3 T — V3 phase diagram from the MD simulations 
showing the freezing and structural transitions. The upper 
and lower arrows correspond to the high and low temperature 
quenches respectively. 



III. MOLECULAR DYNAMICS SIMULATION OF 
SOLID-SOLID TRANSFORMATIONS 

Our aim in this section will be to study the simplest 
molecular dynamics (MD) simulation of the square to 
rhombus transformation. We would like to use the sim- 
ulation results to construct the complete set of coarse 
grained variables needed in a continuum description of 
the dynamics of solid state tranformations. We carry 
out the MD simulation in the constant NVT ensemble 
using a Nose- Hoover thermostat (N = 12000) [§. 

Our MD simulation is to be thought of as a 'coarse- 
grained' MD simulation, where the effective potential is 
a result of a complicated many-body interaction. One 
part of the interaction is a purely repulsive two-body 
potential ^(r^) = v-xjv\^ where r^ is the distance 
between particles i and j. The two-body interaction 
favours a triangular lattice ground state. In addition, 
triplets of particles interact via a short range three-body 
potential V 3 {r t ,r :j ,r k ) = v z w(r i:j ;,r ife ,r ifc )[sin 2 (4% fe ) + 
sin (4:9jki) + sin 2 (49 kij)] where w(r) is smooth short- 
range function and 9ijk is the bond angle at j between 
particles (ijk). Since V3 is minimised when 9ijk = or 
7r/2, the three-body term favours a square lattice ground 
state. Thus at sufficiently low temperatures, we can in- 
duce a square to triangular lattice transformation by tun- 
ing v$. The phase diagram in the T—V3 plane is exhibited 
in Fig. 3. 

We define elastic variables, coarse-grained over a spa- 



Consider two 'quench' scenarios — a high and low tem- 
perature quench (upper and lower arrows in Fig. 3 respec- 
tively) across the phase boundary. In both cases the solid 
is initially at equilibrium in the square phase. 

The high temperature quench across the phase bound- 
ary, induces a homogeneous nucleation (i.e., strain inho- 
mogeneities created by thermal fluctuations are sufficient 
to induce critical nucleation) and growth of a triangular 
region. The product nucleus grows isotropically with the 
size R ~ t 1 / 2 . A plot of the vacancies/interstitial field 
shows that, at these temperatures they diffuse fast to 
their equilibrium value (vacancy diffusion obeys an Ar- 
renhius form D v = Dq exp(— A/ksT), where A is an ac- 
tivation energy and so is larger at higher temperatures). 
The final morphology is a polycrystalline triangular solid. 

The low temperature quench on the other hand, needs 
defects (either vacancies or dislocations) to seed nucle- 
ation in an appreciable time. This heterogeneous nu- 
cleation initiates an embryo of triangular phase, which 
grows anisotropically along specific directions (Fig. 4). 
Two aspects are immediately apparent, the growing nu- 
cleus is twinned and the front velocities are high. In- 
deed the velocity of the front is a constant and roughly 
half the velocity of longitudinal sound. A plot of the va- 
cancy/interstitial field shows a high concentration at the 
parent-product interface. The vacancy field now diffuses 
very slowly and so appear to get stuck to the interface 
over the time scale of the simulation. If we force the 
vacancies and interstitials to annihilate each other, then 



the anisotropic twinned nucleus changes in the course of 
time to an isotropic untwinned one ! 

Therefore the lessons from the MD simulation are : 
(1) There are at least two scenarios of nucleation of a 
product in a parent depending on the temperature of 
quench. The product grows via homogeneous nucleation 
at high T, and via heterogeneous nucleation at low T (2) 
The complete set of slow variables necessary to describe 
the nucleation of solid-solid transformations should in- 
clude the strain tensor and defects (vacancies and dislo- 
cations) which are generated at the parent-product inter- 
face at the onset of nucleation (3) The relaxation times 
of these defects dictate the final morphology. At high 
temperatures the defects relax fast and the grains grow 
isotropically with a diffusive front. The final morphol- 
ogy is a polycrystalline triangular solid. At low temper- 
atures the interfacial defects (vacancies) created by the 
nucleating grain relax slowly and get stuck at the parent- 
product interface. The grains grow anisotropically along 
specific directions. The critical nucleus is twinned and 
the front grows ballistically (with a velocity comparable 
to the sound speed). The final morphology is a twinned 
martensite. 



IV. MODE COUPLING THEORY OF 
SOLID-SOLID TRANSFORMATIONS 

Armed with the lessons from the MD simulation, let us 
now construct a continuum elastic theory of solid-state 
nucleation. The analysis follows in part the theories of 
Krumhansl et. al. H, but has important points of de- 
parture. The procedure is to define a coarse grained free 
energy functional in terms of all the relevant 'slow' vari- 
ables. From the simulation results, we found that every 
configuration is described in terms of the local (nonsin- 
gular) strain field ey, the vacancy field <p, and singular 
defect fields like the dislocation density fry. These vari- 
ables are essentially related to the phase and amplitudes 
of the density wave describing the solid {pa}- 

It is clear from the simulation that the strain tensor, 
defined with respect to the ideal parent, gets to be of 
O(l) in the interfacial region between the parent and the 
product. Thus we need to use the full nonlinear strain 
tensor €y = (diUj+djUi+diiikdjUk)/2- Further, since the 
strain is inhomogeneous during the nucleation process, 
the free energy functional should have derivatives of the 
strain tensor dk£ij (this has unfortunately been termed 
'nonlocal strain' by some authors). 

In general, the form of the free energy functional can 
be very complicated, but in the context of the square-to- 
rhombus transformation, the free energy density may be 
approximated by a simple form, 

f = c(\7e) 2 + 6 2 -ae 4 + e e + Xv4> 2 +Xdb 2 + k d be (1) 

where e is the nonzero component of the strain corre- 
sponding to the transformation between a square and 



a rhombus, <f> is the vacancy field and b is the dislo- 
cation density (we have dropped the tensor indices for 
convenience). The tuning parameter a induces a transi- 
tion from a square (described by e = 0) to a rhombus 
(e = ±e ). 

Starting with e — corresponding to the equilibrium 
square parent phase at a temperature T > T c , we quench 
across the structural transition. The initial configura- 
tion of e is now metastable at this lower temperature, 
and would decay towards the true equilibrium configura- 
tion by nucleating a small 'droplet' of the product. As 
we saw in the last section, as soon as a droplet of the 
product appears embedded in the parent matrix, atomic 
mismatch at the parent-product interface gives rise to 
interfacial defects like vacancies and dislocations. 

Let us confine ourselves to solids for which the energy 
cost of producing dislocations is prohibitively large. This 
would imply that the interfacial defects consist of only va- 
cancies and intcrstitials. The dynamics of nucleation now 
written in terms of e, g (the conserved momentum den- 
sity) and vacancy <\> are complicated B . For the present 
purpose, all we need to realise is that <fi couples to the 
strain and is diffusive with a diffusion coefficient D v de- 
pending on temperature. 

As in the MD simulation, we find that the morphol- 
ogy and growth of the droplet of the product depends 
critically on the diffusion of these vacancies. If the tem- 
perature of quench is high, (j) diffuses to zero before the 
critical nucleus size is attained and the nucleus eventu- 
ally grows into an equilibrium (or polycrystalline) trian- 
gular solid. In this case, the nucleus grows isotropically 
with R ~ t 1 ' 2 . However a quench to lower tempera- 
tures results in a low vacancy diffusion coefficient. In 
the limit D v — > 0, the 0-field remains frozen at the mov- 
ing parent-product interface. In this case a constrained 
variational calculation of the morphology of the nucleus, 
shows that it is energetically favourable to form a twinned 
martensite rather than a uniform triangular structure. 
The growth of the twinned nucleus is not isotropic, but 
along habit planes. Lastly the growth along the longer 
direction is ballistic with a velocity proportional to y/Xv 
(of the order of the sound velocity). All these results are 
consistent with the results of the previous section and 
with martensite phenomenology. Let us try and under- 
stand in more physical terms, why the growing nucleus 
might want to form twins. 

As soon as a droplet of the triangular phase of dimen- 
sion L is nucleated, it creates vacancies at the parent- 
product interface. The free energy of such an inclusion 
is F = Fbulk + Fpp + F(f>- The first term is simply the 
bulk free energy gain equal to AFL 2 where AF is the 
free energy difference between the square and triangular 
phases. The next two terms are interfacial terms. F pp is 
the elastic contribution to the parent-product interface 
coming from the gradient terms in the free energy den- 
sity Eq. [j], and is equal to 4<r pp L, where a pp is the surface 
tension at the parent-product interface. F^, is the con- 
tribution from the interfacial vacancy field glued to the 



parent-product interface and is proportional to <fi 2 ~ L 2 
(since the atomic mismatch should scale with the amount 
of parent-product interface) . This last contribution dom- 
inates at large L setting a prohibitive price to the growth 
of the triangular nucleus. The solid gets around this by 
nucleating a twin with a strain opposite to the one ini- 
tially nucleated, thereby reducing 0. Indeed for an equal 
size twin, <j> — > on the average, and leads to a much 
lower interfacial energy F$ ~ L. However the solid now 
pays the price of having created an additional twin inter- 
face whose energy cost is F tw = a tw L. 

Considering now an (in general) anisotropic inclusion 
of length L, width W consisting of N twins, the free 
energy calculation goes as 



F = AFLW + a pp (L + W) + a tw NW + P[ -^ ) N 



N 



(2) 



where the last term is the vacancy contribution. Mini- 
mization with respect to N gives L/N ~ W 1 ^ 2 , a relation 
that is known for 2-dimensional martensites like In-Tl. 

Our next task is to solve the coupled dynamical equa- 
tions with appropriate initial conditions numerically, to 
obtain the full morphology phase diagram as a function 
of the type of structural change, the parameters entering 
the free energy functional and kinetic parameters like D v . 

It should be mentioned that our theory takes off from 
the theories of Krumhansl et. al. O, in that we write 
the elastic energy in terms of the nonlinear strain tensor 
and its derivatives. In addition we have shown that the 
process of creating a solid nucleus in a parent generates 
interfacial defects which evolve in time. The importance 
of defects has been stressed by a few metallurgists [jj|. 
We note also that the parent-product interface is studded 
with an array of vacancies with a separation equal to the 
twin size. This implies that the strain decays exponen- 
tially from interface over a distance of order L/N. This 
has been called 'fringing field' in Ref. |g{. Krumhansl 
et. al. obtain this by imposing boundary conditions on 
the parent-product interface, whereas here it appears dy- 
namically. 



V. PATTERNING IN SOLID-SOLID 
TRANFORMATIONS : GROWTH AND ARREST 

So far we have discussed the nucleation and growth 
of single grains. This description is clearly valid at very 
early times, for as time progresses the grains grow to a 
size of approximately 1/xm and start colliding, whereupon 
in most alloys they stop. Optical micrographs of acicular 
martensites reveal that the jammed plates lie along habit 
planes that criss-cross and partition the surrounding fee 
(parent) matrix. 

Can we quantify the patterning seen in martensite ag- 
gregates over a scale of a millimeter ? A useful measure 



is the size distribution of the martensite grains embedded 
in a given volume of the parent. The appropriate (but 
difficult !) calculation at this stage would be the analogue 
of a Becker-Doring theory for nucleation in solids. In the 
absence of such a theory, we shall take a phenomcnolog- 
ical approach. 

Clearly the size distribution P(l,t) depends on the 
spatio-temporal distribution I of nucleation sites and the 
growth velocity v. We have analysed the problem explic- 
itly in a simple 2-dimensional context. Since the nucle- 
ating martensitic grains are highly anisotropic and grow 
along certain directions with a uniform velocity, a good 
approximation is to treat the grains as lines or rays. 
These rays (lines) emanate from nucleation sites along 
certain directions, and grow with a constant velocity v. 
The rays stop on meeting other rays and eventually after 
a time T, the 2-dimensional space is fragmented by N 
colliding rays. The size distribution of rays, expressed in 
terms of a scaling variable y — y(I,v), has two geomet- 
rical limits — T-fixed point (at y — 0) and the L-fixed 
point (at y = oo). The T-fixed point corresponds to 
the limit where the rays nucleate simultaneously with a 
uniform spatial distribution. The stationary distribution 
P(l) is a Gamma distribution with an exponentially de- 
caying tail. The L-fixed point, corresponds to the limit 
where the rays arc nucleated sequentially in time (and 
uniformly in space) and grows with infinite velocity. By 
a mapping onto a multifragmentation problem, Ben Nairn 
and Krapivsky p0[ were able to derive the exact asymp- 
totic form for the moments of P(l) at the L-fixed point. 
The distribution function P(l) has a multiscaling form, 
characterised by its moments < l q >~ N-v(i) where 
/■*(<?) = (? + 2 — \Jq 2 + 4)/2. At intermediate values of 
the scaling variable y, there is a smooth crossover from 
the T-fixed point to the L-fixed point with a kinematical 
crossover function and crossover exponents. 

The emergence of scale invariant microstructures in 
martensites as arising out of a competition between the 
nucleation rate and growth is a novel feature well worth 
experimental investigation. There have been similar sug- 
gestion in the literature, but as far as we know there has 
been no direct visualization studies of the microstructure 
of acicular martensites using optical micrographs. Re- 
cent acoustic emission experiments [ pT[ on the thermoe- 
lastic reversible martensite Cu-Zn-Al, may be argued to 
provide indirect support of the above claim , but the 
theory of acoustic emission in martensites is not under- 
stood well enough to make such an assertion with any 
confidence. 



VI. OPEN QUESTIONS 

We hope this short review makes clear how far we are 
in our understanding of the dynamics of solid-solid trans- 
formations. A deeper understanding of the field will only 
come about with systematic experiments on carefully se- 



lected systems. For instance, a crucial feature of our 
nonequilibrium theory of martensitic transformations is 
the existence of a dynamical interfacial defect field. In 
conventional Fe based alloys, the martensitic front grows 
incredibly fast, making it difficult to test this using in 
situ transmission electron microscopy. Colloidal solutions 
of polysterene spheres (polyballs) however, are excellent 
systems for studying materials properties. Polyballs ex- 
hibiting fee — > bee structural transitions have been seen 
to undergo twinned martensitic transformations. The 
length and time scales associated with colloids are large, 
making it comfortable to study these systems using light 
scattering and optical microscopy. 

In this article we have focussed on a small part of the 
dynamics of solid state transformations, namely the dy- 
namics and morphology of martensites. Even so our pre- 
sentation here is far from complete and there are crucial 
unresolved questions that we need to address. 

Let us list the issues as they appear following a nucle- 
ation event. 

The physics of heterogeneous nucleation in solids is 
very poorly understood. For instance, it appears from 
our simulations that the morphology of the growing nu- 
cleus depends on the nature of the defects seeding the 
nucleation process (e.g., vacancies, dislocations and grain 
boundaries). In addition several martensitic transforma- 
tions are associated with correlated nucleation events and 
autocatalysis. Though these features are not central to 
the issue of martensites, such a study would lead to a bet- 
ter understanding of the origins of athermal, isothermal 
and burst nucleation. This in conjunction with a 'Becker- 
Doring theory' for multiple martensite grains would be a 
first step towards a computation of the TTT curves. 

We still do not understand the details of the dynamics 
of twinning and how subsequent twins are added to the 
growing nucleus. Moreover the structure and dynamics of 
the parent-product interface and of the defects embedded 
in it have not been clearly analysed. 

It would be desirable to have a more complete theory 
which displays a morphology phase diagram (for a single 
nucleus) given the type of structural transition and the 
kinetic, thermal and elastic parameters. 

Certain new directions immediately suggest them- 
selves. For instance, the role of carbon in interstitial 
alloys like Fe-C leading to the formation of bainites; the 
coupling of the strain field to an external stress and the 
associated shape memory effect ; and finally the nature of 
tweed phases and pre-martensitic phenomena (associated 
with the presence of quenched impurities). 

It is clear that the study of the dynamics of solid-solid 
transformations and the resulting long-lived morpholo- 
gies lies at the intersection of metallurgy, material science 
and nonequilibrium statistical mechanics. The diversity 
and richness of phenomena make this an extremely chal- 
lenging area of nonequilibrium physics. 
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Figure Caption 

Fig. 4 MD snapshot of (a) the nucleating grain at 
some intermediate time initiated by the low tempera- 
ture quench across the square-triangle transition. The 
dark(whitc) region is the triangular (square) phase re- 
spectively. Notice that the nucleus is twinned and highly 
anisotropic, (b) the vacancy (white) /interstitial (black) 
density profile at the same time as (a). Notice that the 
vacancies and interstitials are well separated and cluster 
around the parent-product interface. 
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